Paired composite fermion phase of quantum Hall bilayers at 
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We provide numerical evidence for composite fermion pairing in quantum Hall bilayer systems at 
filling ^ = I + 5 for intermediate spacing between the layers. We identify the phase as p^ + ipy 
pairing, and construct high accuracy trial wavefunctions to describe the groundstate on the sphere. 
For large distances between the layers, and for finite systems, a competing "Hund's rule" state, or 
composite fermion liquid, prevails for certain system sizes. 



The bilayer quantum Hall system at filling fraction 
V = ^ +_i has been an active topic of research for over a 
decade [ll Q . In the limit of small layer spacing d there 
is a superfluid phase, whereas in the limit of infinite d 
the two layers form independent compressible composite 
fermion (CF) seas H- Much less well understood, how- 
ever, is the nature of the state at intermediate d. Even in 
simplified models with spinless electrons at zero temper- 
ature, no disorder, and no tunneling between layers (sim- 
plifications we will adopt throughout this paper), there 
are a wide range of theoretical calculations, predictions, 
and scenerios [J, la Is, LZ| attempting to address the transi- 
tion. One of the most interesting possibilities [6| is that, 
even for very weak interaction between layers, the inter- 
action between the two CF seas may cause BCS pairing 
between the layers to form a superconducting state: an 
"interlayer paired CF-BCS state" . Unfortunately, this 
suggestion was based on a highly approximate pertur- 
bative Chern-Simons approach whose validity at higher 
order is impossible to test. In the current paper we pro- 
vide highly accurate numerical evidence that for a range 
of distances between the layers, the groundstate is in- 
deed such a CF-BCS state. We further show that the 
pairing symmetry is in the Px + ipy channel in contrast 
to the theoretical arguments of Ref. |5| which suggest 
Px — iPy We construct and numerically verify the first 
explicit wavefunctions describing this phase. For finite 
systems for very large d we also identify the groundstate 
to be a "Hund's rule" or composite fermion hquid state. 
For certain shell fillings (number of particles) this is pre- 
cisely the weak-interaction limit of the CF-BCS state, 
whereas for other fillings, it is a distinct state. Finally, 
we argue that, when the two possibilities are distinct, the 
range in d for which the Hund's rule state prevails over 
the CF-BCS state becomes smaller as we go to larger sys- 
tems while the paired state becomes more predominant. 

Since many of our technical details are similar to that 
of Ref. ja|, where single-layer pairing is considered, we 
will be brief with our discussion. We start by consider- 
ing the wavefunction J|^ ( 1 + <7k ^^^ '^kT -ki ) l'^) which 
represents a BCS paired wavefunction [9| for a bilayer 
spinless fermion system in zero magnetic field (where 



t , J, represent the layer index) . Fourier transforming 
with respect to Lp we project to precisely N /2 fermions 
(TV assumed even) in each layer [9[, yielding ^bcs — 
det [G(rj,| - Yj^i)] where G(r,,| - r^j) is an N/2 by N/2 
matrix with indices i and j. G(ri^| — Vj^i), the wavefunc- 
tion of a pair, can be written in terms of gt as 



G(r»,T -rj,i) ^ Ek.9k0k(rj,T)'/'-k(rj,i) 



(1) 



where (/)k(r) — e^^'^ are the simple single fermion plane 
wave orbitals, such that Eq. [T] is just a Fourier trans- 
form. BCS theory is fundamentally variational [9|, and 
one uses a form of ^k which minimizes the energy given a 
particular inter-particle interaction. The "symmetry" or 
relative angular momentum of the pairing wavefunction 
is determined by the phase winding [10] of G. If the phase 
of G(ri.| — Tj^i) wraps by 2M-K as r^^j is taken clockwise 
around Yj^i, we say the wavefunction has angular mo- 
mentum M (or M-wave symmetry) where M = 0, ±1 are 
also known as s,px^ ipy respectively. (While some prior 
literature uses "pa; -I- ipy^ imprecisely to denote either 
chirality). Note that the limit of noninteracting fermions 
can be achieved with this form by taking (?k > for all 
|k| < kp and gt = otherwise (with kp the Fermi mo- 
mentum). Thus, the BCS paired state can be deformed 
smoothly into a noninteracting Fermi liquid (a limit point 
outside of the superconducting phase). 

Our calculations are performed on a spherical ge- 
ometry. We remind the reader that with a magnetic 
monopole of flux 2q at the center of the sphere (with 
q half integer), the single particle eigenstates [ll| are 



the spherical monopole harmonics Yl^^{r) where I — 
q,q + \, . . ., and m = —I, —l-\-l,...^l and the energy 
of these states depend on / only. From these orbitals we 
construct a general pair wavefunction 



G(r 



i^T-rj 



i) = E,™5K-i)'+">^:'™(r,T)r,%„(r,a) (2) 



where again the gi are variational parameters that con- 
trol the (radial) shape of the pair wavefunction. It is 
easy to show that this pairing wavefunction has M-wave 
symmetry with M = 2q. Indeed, one can show that in 
order to have a BCS groundstate with M-wave symme- 
try on the sphere with no vortex defects, one must have 



a monopole of flux M at the center of the sphere, and the 
pairing form is always expressible in the form of Eq. [2] 

To convert our siinple BCS wavefunction into a paired 
composite fermion [3|, Il2| wavefunction appropriate for 
modeling the bilayer quantum Hall state at v = -^ -\- -^^ 
we must multiply by Jastrow factors such that electrons 
"see" flux attached to electrons only within the same 
layer. To achieve this, we replace the single particle or- 
^ with composite fermion orbitals defined as 

%,a,Vi,a) = 7^[Ji,ai^',„(i"j,a)] whcrc Wi,Q,Wi,Q is the 
usual spinor representation of the coordinate Vi^a on the 
sphere, a =t, i is the layer index, Ji^a = Y\j^i{ui,aVj^a — 
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,) is the composite fermionization factor, and V is TABLE I: Data for bilayer electrons on a sphere near v = \ + l 



the projection to the lowest Landau level. Note that in 
our present notation, the Jastrow factor J is absorbed 
within the composite fermion wavefunction Y (as op- 
posed to Refs. H^])- As usual, Yi\^{ui^a,Vi^a) is imphc- 
itly a function of all the particle coordinates, although 
we only denote it explicitly as a function of particle i. 

We substitute these composite fermion orbitals Y in 
place of the orbitals Y in Eq. [2] to generate the pair- 
ing wavefunction which we will correspondingly call G. 
Our trial wavefunction for the bilayer quantum Hall 



state is then VP = det 



completely analo- 



gous to the case of simple BCS theory discussed above. 
The generated trial wavefunction (for N/2 electrons in 
each layer) occurs on the sphere with a monopole flux 
N^ = 2{N/2 — 1) + M for the case of M-wave symmetry 
(this "shift" of AI is caused by the small addition of M 
flux quanta necessary to avoid having vortices, as dis- 
cussed above) . The different possible pairing symmetries 
can then be easily identified by their shifts. 

Our numerical analysis is founded on exact diago- 
nalizations on the sphere. The inter-electron interac- 
tion is taken to be V{ri^a,^j.p) = e^/r ioi a — (3 and 
e^/Vr"^ + <P for a 7^ /? where d represents the "spacing" 
between layers (a,/3 =t,i are the layer indices). We 
have defined the magnetic length to be unity, and r is 
the chord-distance between r^ and Yj . In table U we ex- 
amine the stability of the bilayer system at different shifts 
(different values of the magnetic monopole flux) . We re- 
call that quantum Hall states correspond to groundstates 
at angular momentum L = which have a strong gap. 
While it is obvious that there should be a gap and an 
L = groundstate when each individual layer is gapped 
with an L = groundstate (which is the case for A^ = 12 
for both M = +1,-1), it is less trivial in the other cases 
where the individual layers have L ^ groundstates. 
We find that the M = +1 case has an L = ground- 
state with a strong gap for all values of the interlayer 
spacing d, while this is not the case for other values of 
M {M = 0,±1,±2 have been examined). We thus sus- 
pect that if a BCS paired state of CF's exists, it is in the 
Px + ipy channel (M = +1). We will thus focus on trial 
CF-BCS wavefunctions with this pairing symmetry. 



for various interlayer spacings d (in units of the magnetic 
length £o). The table shows energy gaps (in units of e'^/^o) for 
cases where the groundstate is angular momentum zero (non- 
bold), and angular momentum L of the groundstate (bold) 
when it is not zero. Data is shown with N/2 electrons per 
layer with flux Ai'^ = iV - 2 -I- M for A/ = -1, 0, -f 1 corre- 
sponding to the pairing symmetries px — ipy, s, and p^ + ipy. 
Only the case of px + ipy consistently shows a strong gap 
at zero angular momentum characteristic of a quantum Hall 
state. Data for M = ±2 (not shown) also does not suggest 
a quantum HaU state. We thus identify Px + ipy as the most 
likely pairing symmetry. Note that the strong gap for A'^ = 12 
in the M = —1,-1-1 case both correspond to flUed shells of 
composite fermions. In non-filled shells, the gap drops with 
increasing d as the interlayer interaction is reduced. 



Using the above described approach, we generate trial 
CF-BCS wavefunctions with p^ + ipy pairing. Calcula- 
tions are performed by Monte-Carlo, and at each value of 
the interlayer spacing d, the shape of the pairing wave- 
function is optimized by varying the parameters gi to 
maximize the overlap with the exact groundstate (See 
Ref. [8| for details of the calculational scheme). Very 
good agreement with the exact diagonalization is ob- 
tained by only varying the first few such parameters. 
In Fig. [T^ we show overlaps of our trial wavefunctions 
with the exact groundstate for several different sized sys- 
tems {N = 10, 12, 14) over a range of d. For d > 1 the 
overlaps of the trial state with the exact groundstate are 
excellent. Since these wavefunctions are variational, one 
might worry that a large number of variational param- 
eters is required to obtain such good agreement. How- 
ever, the dimension D of the L ^ Hilbert space is 
D = 38, 252, 1559 for N = 10, 12, 14, and we have used 
only 3, 4 and 4 variational parameters gi respectively, so 
these agreements are highly significant. 

Also shown in Fig. [T^ is the overlap of the exact 
groundstate with the CF-Fermi sea trial wavefunction, 
as well as that with the Ill-state. For the iV = 12 case, 
the CF-Fermi sea state is uniquely defined as filling two 
CF shells in each layer (at this flux, the p*'' CF shell 
has 2p orbitals) [3|, \l^ . At large d this CF-Fermi sea has 
virtually perfect overlap with the exact groundstate. Fur- 
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FIG. 2: (color online) Energy difference between CF-BCS 
wavefunction and Hund's rule wavefunction as a function of 
layer spacing d for shell fillings where CF-BCS and Hund's 
rule are distinct states in the large d limit. We show A^ = 
8, 16, 20. Note that as A'^ grows, the range of d where the 
CF-BCS wavefunction is a better trial state extends to larger 
d. The inset shows dc, the value of d where gbcs ~ enund to 
illustrates how it scales with N~^; lines are guides to the eye. 



FIG. 1: (color online) Overlaps of our trial CF-BCS states 
with the exact groundstate on the bilayer sphere (solid), over- 
laps of the Hund's rule or CF-Fermi sea state with the exact 
groundstate (dashed) and overlaps with the Ill-state (dot- 
ted). Error bars are Monte-Carlo error. In (a) data is shown 
for A'^ = 10, 12, 14 where the large d limit of the Hund's rule 
and CF-BCS state become identical. Note the very high ac- 
curacy of these trial wavefunctions for all d > 1. In (b) for 
A^ = 8, 16 the two limits are inequivalent. In this case it is 
clear that the Hund's rule state becomes more accurate at 
large d, but at smaller d there is a regime where the CF-BCS 
state prevails. The Ill-state is accurate only at d < 0.5^o. 



thermore, the bilayer CF-Fermi sea in this case is simply 
a hmit of the CF-BCS wavefunction (as described above). 
However, at smaller d, where the CF-BCS wavefunction 
achieves almost perfect agreement with the exact ground- 
state, the CF-Fermi sea has very poor overlap. 

For values of N which are not filled shell situations, 
the CF-Fermi sea needs to be more carefully defined. We 
recall that for a single-layer system (at least for small sys- 
tems) , the groundstate of a partially filled shell satisfies 
Hund's rule by filling degenerate CF-orbitals to maxi- 
mize the total angular momentum [13[ . We thus propose 
that the groundstate of the bilayer at very large d is ob- 
tained by combining two single- layer states, each of which 
satisfies Hund's rule individually, into an overall angular 
momentum singlet (adding up the wavefunctions of the 
two individual layers into the unique Ltotai = state 
using appropriate Clebsch-Gordan coefficients). We call 
this construction the "Hund's rule" bilayer state. We 
see in Figs.[T^ and[T)D that the Hund's rule state always 
becomes extremely accurate at large d. In cases where 



there is either one CF in the valence shell in each layer 
(such as iV = 14) or one CF-hole in the valence shell in 
each layer(such as iV = 10), it can be shown that the 
Hund's rule state is again precisely the large d limit of 
the CF-BCS state. Indeed, the Hund's rule and CF-BCS 
state become identical in this limit (and both become al- 
most perfect). However, in other cases where there is an 
incompletely filled shell {N = 8,16 shown in Fig. [T]]), 
the Hund's rule state cannot be written as a limit of the 
CF-BCS state. To understand that these must be differ- 
ent, we note that for the Hund's rule state, each layer 
is an eigenstate of L^, whereas for the BCS state this is 
only true in the limit where the pairing becomes infinitely 
weak and when it is a filled shell situation, or there is a 
single electron or single hole in each valence shell. 

As seen in Fig. [T]d, for A^ = 8, 16, the Hund's rule 
state becomes extremely accurate in the large d limit, 
and the CF-BCS state becomes inaccurate. Remarkably, 
as we go to smaller d, the CF-BCS state again becomes 
extremely accurate whereas the Hund's rule state fails. 
This is an extremely important result: even when the 
large d limit is not of the CF-BCS pairing form, when the 
interaction between layers is increased, the pairing form 
again becomes accurate. (For iV = 8, 16 the dimension 
of the i = Hilbert space is D — 12, 1277 A and we have 
used 2,4 variational parameters, so again the agreement 
is significant. In all cases, the Ill-state which describes 
the interlayer coherent phase is accurate only at small 
d < 0.5^0 J where its overlap with the exact groundstate 
plummets and the CF-BCS states become the best trial 
states. Thus, we argue that there is a region of interme- 
diate d where a CF-BCS phase is the groundstate. 

We further conjecture that as we go to larger systems, 



the shell-filling effects, and Hund's rule, should become 
less important whereas the pairing effects will remain 
the same strength. This conjecture is simply based on 
the fact that only ~ ^/N particles are in the valence 
shell, whereas all particles within some gap energy of 
the Fermi surface (a number ^ N) contribute to pair- 
ing. To make this statement more concrete we examine 
the range for which the CF-BCS wavefunction provides a 
better trial state than the Hund's rule state. In Fig. [2] we 
show the energy difference per particle between the two 
trial states as a function of system size and layer separa- 
tion. To differentiate between the two possibilities most 
clearly, we have only shown cases {N = 8, 16, 20) where 
the large d limit of the CF-BCS state is distinct from the 
Hund's rule state. The variational CF-BCS wavefunc- 
tion for A^ = 20 electrons was obtained using an energy 
minimization technique 1J| . The figure shows that with 
increasing N , the CF-BCS state becomes more accurate 
out to somewhat larger d. While our data strongly sug- 
gests that pairing survives in the thermodynamic limit it 
cannot establish that it extends to arbitrary weak inter- 
actions, as suggested in Ref. [6|, in this limit (see inset). 

For large d, it is difficult to be confident of extrapo- 
lation to the thermodynamic limit (small energy differ- 
ences). However, for intermediate io ^ d <^ 2£o, our 
data (table |T]) clearly shows a strong gap which persists 
to large N . At these values of d, it is also clear from 
our numerics that the CF-BCS state provides an accu- 
rate trial wavefunction which is not in the 111 phase. We 
can thus conclude that a CF-BCS phase, contiguous to 
the 111 phase, does exist for a range of intermediate d. 

We emphasize that the CF-BCS trial wavefunctions 
contain interlayer correlations via the interlayer pairing 
function G, whereas the intra-layer correlations result 
from the CF Jastrow factors. As such, these trial wave- 
functions should not have the same interlayer coherence 
as the 111 wavefunction. Still, our CF-BCS states should 
have quantized Hall drag exactly like the 111 phase, at 
least at zero temperature (see [5|, Il5[), although they 
should not display resonant interlayer tunneling. Finally, 
we note that, unlike the 111 phase, the CF-BCS state will 
be strongly destabilized by layer imbalance (analogous to 
spin imbalance for usual BCS states). 

In a forthcoming paper 14| we will discuss the overall 
phase diagram including the regime of very small d. In 
brief, we find a second order phase transition near (i > 1 
from the CF-BCS phase at larger d, which has zero 111 
order parameter (OP), to a phase with nonzero 111 OP 
at smaller d. Note that the BCS OP and the 111 OP 
are distinct. We find that both phases as well as the 
transition can be very well described in the language of 
[J| where CF's mix with composite bosons (CB's) and the 
presence of CB's yields nonzero 111 OP. Our work further 
suggests that there may be a region of intermediate d 
where both OPs (BCS and 111) coexist, although the 
results of 14| are not definitive in this respect. Should 



this coexistence occur, we argue [ij] that px -\- ipy is the 
only pairing channel compatible with such coexistence. 

To summarize, we have shown compelling numerical 
evidence of px + ipy pairing of composite fermions at 
intermediate layer spacings d for quantum Hall bilayers 
at i^ = 5 + 5- We have proposed specific forms for the 
actual wavefunctions that show excellent overlap with 
the results of exact diagonalizations. While CF pairing 
had been theoretically proposed earlier y, other phases 
had also been advocated [5|, LZf and there has previously 
been no compelling numerical evidence to distinguish the 
possibilities. Further we show that for finite size systems 
at very large d, a Hund's rule state is the groundstate. 
However, as we go to larger system sizes, the CF-BCS 
(paired) state extends to larger d. 
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